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THE FRONT LOCATION IN BBM WITH DECAY OF MASS 


LOUIGI ADDARIO-BERRY AND SARAH PENINGTON 


Abstract. We augment standard branching Brownian motion by adding a competitive 
interaction between nearby particles. Informally, when particles are in competition, the 
local resources are insufficient to cover the energetic cost of motion, so the particles’ 
masses decay. In standard BBM, we may define the front displacement at time t as the 
greatest distance of a particle from the origin. For the model with masses, it makes sense 
to instead define the front displacement as the distance at which the local mass density 
drops from 0(1) to o(l). We show that one can find arbitrarily large times t for which 
this occurs at a distance Oft 1 / 3 ) behind the front displacement for standard BBM. 


1. Introduction 

In this work, we propose a mathematical model of competition for resources within a 
single species, in a growing, spatially structured population, and provide an initial study of 
the front location in this new setting. The model is essentially standard one-dimensional 
branching Brownian motion (BBM), augmented with a destructive, local interaction be¬ 
tween particles. We first briefly recall BBM: start from a single particle at a point in R, 
endowed with an Exp(l) “branching clock”. The particle moves according to Brownian 
motion; when its clock rings, it splits in two (branches). The new particles receive inde¬ 
pendent Exp(l) clocks, and move independently (according to Brownian motion) starting 
from where the first particle splits, until their own clocks ring and they in turn split, et 
cetera. 

Write n(t) for the total number of particles at time t, and X(t) = (Xj(t), 1 < i < n(t)) 
for the locations of such particles. We assume the particles are listed in a way that makes 
the vector X(t) exchangeable; one possible formalism is via the Ulam-Harris tree, with 
particles listed lexicographically acording to their label in the tree. We refer the reader to 
[10] for more details on such matters; but many different references are possible. We also 
write N(t,x ) = {i : Xi(t ) > x\ for the indices of particles with position greater than x at 
time t. 

We sometimes write > 1), ignoring the fact that X(t) has finite length, for 

convenience. We adopt the convention that Xp(t) = d for k > n(t) (so d is where new 
babies come from). We refer to “the particle X y(t)” as shorthand for “the particle with 
position Xi(t) at time t”; this is unambiguous at Leb-a.e. time t. We write P for the 
probability measure under which (X(f),f > 0) has the law of one dimensional BBM with 
initial individual at 0, E for the corresponding expectation, and (J-), t > 0) for the filtration 
generated by the process. 

We now add destructive interaction as follows. Informally, imagine that the particles 
are, say, amoeba. Motion has an energetic cost, but for a single particle in isolation, this 
cost is exactly accounted for by the resources available in the environment. When particles 
are nearby, however - at distance less than one, say - they must share resources; in this 
case individuals do not consume enough to meet their energy expenditure, and their mass 
decreases. Finally, larger (more massive) individuals consume resources at a greater rate. 
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Formally, we define a vector M(f) = (M*(f),i > 0), and call M % (t) the mass of particle 
Xi{t). By convention, if X*(f) = d then M*(f) = 0. Write 

C (t,x)= y M *( s ) 

{i:|Xi(t)-x|e(0,l)} 

for the total mass of particles within distance one of x at time t, excluding any particles 
at position x. Then at time t, Mi(t ) decays at rate £( t,Xi(t )). In other words, dMj(t) = 
-Mi(t) ■ C(t,Xi(t))dt, so 

Mi(t) = exp(— [ C(s,Xj(s))ds). 

Jo 

This should be viewed as defining (M(f), t > 0) to be the solution of a system of differential 
equations; the definition makes sense since the system has a unique solution P-almost 
surely. Furthermore, the process (M(f),f > 0) is clearly J^-adapted. 

For later use, write (Xj^s),!) < s < t) for the ancestral path leading to X*(f), and let 
Oi{t) be the final branching time along this path. Also, let ji,t(s) be the index of Xj it (s) 
among the time-s population, so that Xj )i (s) = Xj. t ( s )(s). We also write Mij(s) for the 
mass of the ancestor of X*(f) at time s (so Mi :t (s ) = Mj. t ( s )(s)). 

Note that along any given trajectory, mass decreases: (M^(s),0 < s < t) is decreasing 
in s for each 1 < i < n(t). Mass enters the system through branching events, since each 
“child particle” inherits the mass of its parent. This is obviously physically unrealistic 
in some settings (e.g. for amoebae) but may be more realistic in others (e.g. in nuclear 
physics). 

Rather than viewing Afj(t) as a mass, a perspective suggested by a referee is to view 
(Mj^(s),0 < s < t) as recording information about the local density of the environment 
observed along the ancestral trajectory of the particle Xj(t). The interaction between the 
dynamics of X(f) and M(f) makes this point of view slightly complicated to interpret, 
but here is one possibility. Imagine adding destructive interaction to a BBM, as follows: 
whenever two different particles are at distance less than 1, each kills the other at rate one. 
Record such a killing event as a mark at the appropriate location of the BBM family tree. 
Particles with a mark on their ancestral trajectory are ghosts, which continue to move and 
reproduce as before, but can no longer kill other particles. Given the BBM but not the 
marks, one may ask for the conditional survival probabilities Pi(t) = P {Xj(t) is alive | J~t } 
of the particles. The vector M(t) is a “linearized” version of the vector of these survival 
probabilities. 


l. 1. Main result. Write 

d{t, m) = minjx > 0 : £(f, x) < m}, D(t, m) = maxjx : ((t , x) > m}, 

for the leftmost (positive) location at which the total mass of nearby particles falls below 

m, and the rightmost location at which it exceeds m, respectively. We prove the following 
theorem. 


Theorem 1.1. Write c* = 3 4 / 3 tt 2 / 3 /2 7 /6 . Then almost surely, for all m < 1, 

and 


V2 t-d(t,m) * 

inn sup - “tt: - > c 

t-> oo t 1 / 3 


Uminf 

£—>•00 £ 1 /^ 


A well-known result of Bramson [4] states that the rightmost particle location maxj>i Xj(t) 
has median med(t) satisfying 

med(t) = \/2t — logt + 0(1). 

Furthermore, it turns out [12] that | maxj>i Xj(t)— med(f)| is almost surely 0(log t), in that 
limsup^gQ | maxj>i X;(t) — med(f)|/log t is a.s. finite. In view of this, the theorem states 
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that (1) there are (large) times t at which the first low-density region lags at least distance 
c*t 1//3 + o(t 1 / 3 ) behind the rightmost particle, and (2) there are also (potentially different, 
large) times t at which there is some high-density region within distance c*f 1//3 + o(t 1 / 3 ) 
of the rightmost particle. 

We believe that in fact almost surely, for all m < 1, 


lim 

t->-oo 


y/21 — d(t, m) 

3 


= c* = hm 

£—>•00 


\/2 1 — D(t, m) 
tV3 


If this is correct, then the front could equivalently be defined as, e.g., a median of D(t,m ) 
or d(t, m) - or any other fixed quantile of one of these random variables. We provide some 
justification for our belief in Section 6. That section also contains a few open questions 
about the model and a discussion of various generalizations of our results (some straight¬ 
forward, some conjectural), as well as describing variants of the model which have thus 
far resisted analysis. 


2. Proof sketch 

Here comes an outline of the key tools in our argument. The first is technical but impor¬ 
tant and also, we believe, provides important intuition when making heuristic predictions 
about the behaviour of the process. The remainder gives a fairly detailed overview of the 
proof. 

Density self-correction. It is not hard to see that when £(t, x) is small (much less than 
one), and this also holds in a region around x, then £(t, •) will exhibit exponential growth 
near x, at least for a short time. Indeed, we heuristically have 

^C(M) ~ C{t,x) - ^2 Mi(t)-((t,Xi(t)). 

This is not exactly correct since it ignores the effect of motion (particles may enter or leave 
the region near x), but it is a useful first approximation. In particular, it suggests that 
if £(t, y) is small (much less than one) for all y with \y — x\ < 1, then £(f, •) will exhibit 
exponential growth near x, at least for a short time. This is indeed true; one important 
consequence is that if £(i, x) = e and C(V) is not too wild then it is very likely that 
C(t',x) = 0(1) for some t' = t + @(log(l/e)). Similarly, when £(t,y) is much larger than 
1 for y near x then £(t, x) will decrease exponentially quickly. We use the self-correcting 
nature of the density in several places throughout the paper. 

As an aside, we remark that if ((t,y) ~ £( t,x ) for \y — x\ < 1 then the above heuristic 
gives 4fC{t,x) ~ £(t, x)( 1 — C(t,x)), which is suggestive of the logistic control; we briefly 
revisit this connection in the conclusion. 

Population + no competitors=mass. Fix /3 > 0 and suppose that for some function 
/ : [0, oo) — i K, for all s € [0,i], D(s,/3) < f(s), or in other words ((s,x) < j3 for all 
x > f{s). In this case, particles that stay ahead of the moving barrier / are in a relatively 
sparse environment, so do not lose mass too quickly. More precisely, if Xi(t) satisfies 
Xi t t(s) > f(s ) for all s € [0, t] then Mj(f) > e _/3t . It follows that for any x > f(t ) + 1, 

C (t,x) > e~ !3t •#{* : I Xi(t) -x\< 1, Vs € [0, f], X ijt (s) > /(s)}. 

For such x, if #{i : | Xi(t) — x\ < l,Vs € [0, t], Xij(s) > /(s)} > /3e^ then £(i, x) > f3, 
contradicting the assumption that D(t,j3) < f(t). 

Surfing the wave. To exploit the above contradiction, we require that with high prob¬ 
ability there are many particles staying ahead of some barrier. Such results are available: 
it follows fairly straightforwardly from recent studies of consistent maximal displacement 
for BBM [20] that for c > c*, for all large times t there are e®^ 1/3 ) particles at time t 




4 


LOUIGI ADDARIO-BERRY AND SARAH PENINGTON 


with which have stayed ahead of the curve f(s) = \/2s — cs 1 / 3 . This allows us to take 
/3 = t _1 above and obtain that there is s € [0,f] and x > f(s) such that C( s > x ) > t ■ 
Since the local density grows exponentially in regions with small density, we will with 
high probability find s' with £(s',x) > b > 0 and s' — s = O(logt). For such s' we 
have f(s') = f(s ) + O(logt) so x > V2s' — c(s') 1 / 3 — O(logt) so with high probability 
d(s , ,b ) > y/2s' — c(s') 1/3 — O(logf). 

The lower bound is practically complete, but we must rule out the possibility that 
s' = 0(1) for all t. To do so, we first establish that 


sup 
t> o 


ma x{((t,x),x € M} 
log(t + 2) 


= : Z < oo 


almost surely. 


Proving this is harder than might be expected; its proof, given in Section 4, occupies 8 
pages and is perhaps the most technically challenging part of the paper. 

Once we prove that Z < oo, we then reprise the above argument, but with a variable 
mass bound 


p = m 


Z log(s + 2) for s < f 1 / 4 
t~ l for s € (t 1 / 4 , t] . 


The loss of mass before time t 1 / 4 is insignificant compared with that which follows, so 
essentially the same argument as above yields that there is s € [0, t] and x > f(s ) such 
that C(s,x) > P(s). On the other hand, this can not happen for s < t 1 / 4 by the definition 
of Z, so it must happen later. This is enough to conclude the lower bound. The details of 
this argument appear in Section 5. 


Competition implies decay. For the upper bound, given in Section 3 (with some 
technical lemmas deferred to an appendix), we invert the above argument by contradiction. 
In brief: if all particles to the right of a given curve have spent large amounts of time in 
high-mass environments, then all such individuals will have very low weight; if furthermore 
there are not many of them, then their total weight is also small. 

More precisely, suppose that for some C > 0 and some function g : [0, oo) —> M, for all 
s € [CfV 3 /2,f], we have d(s,m ) > g(s), so ( (s,x) > m for all x with x € [0, g(s)]. Then 
for all i, 

Mi(t ) < exp(-m • Leb(s € [Ct 1/3 /2, t] : |X M (s)| € [0,s(s)])). 

It follows that if all particles with Xi(t) > g(t ) have Leb(s € [Ct 1 / 3 /2,t] : |Xj )t (s)| € 
[0, g(s)]) > l then for all x > g(t) +1, recalling the notation N(t, x) from the introduction, 

C(t,x)<e- mi -\N(t,gm- 

If \N(t, <?(t))| < me mi , this is in contradiction with the assumption that d(s,m) > g(s). 


Whitecaps are just foam. Once again using estimates related to consistent maximal 
displacement for BBM, we show that for c < c*, with g(s) = y/2s — cs 1 ^ 3 , with high 
probability every particle with Xi(t) > g(t) indeed spends at least a time Ct 1 ^ 3 behind the 
curve g. Under the above assumption, it follows that the particles counted by N(t,g(t )) 
are as insubstantial as sea spray; for all x > g(t) + 1, 

C (t,x) < e~ mCtl/3 / 2 • #{i : I Xi(t) - x\ < 1} < e~ mCtl/3 / 2 ■ \N(t,g(t))\- 

Standard and simple arguments for BBM show that \N(t, g(t))\ = with high prob¬ 

ability, so we obtain a contradiction for large t if C is sufficiently large. It follows that 
with high probability there is s € [Ct 1 / 3 /2, t] and x € [0,g(s)] such that ((s,x) < m: for 
such s we have d(s,m) < g(s). This is the content of Proposition 3.2. 
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Definitions 

We sometimes need to consider the evolution of a subset of the particles starting at a 
time greater than zero, so it is useful to allow initial conditions other than a single mass- 
one particle at the origin. Generally, for x = (xi,. .., Xf.) G M. k and m = (mi,..., m^) G 
(0, oo) k , we write P x . m for the probability measure corresponding to an initial condition 
with a particle of mass m; at location x* for each 1 < i < k. We write P = P(o),(i) f° r the 
default initial condition. 

We say a random variable X is geometric with parameter p, or is Geom(p)-distributed, 
if P {X = &} = (! — p) k ~ 1 p for positive integer p. 


3. Upper bound 

Recall from the introduction that c* = 3 4 / 3 7 t 2 / 3 2~'/ 6 . The next proposition is a restate¬ 
ment of the upper bound from Theorem 1.1. 


Proposition 3.1. For any m > 0, almost surely 


lim sup 

£—>■00 


y/2t — d(t , m) 
tV3 


> c*. 


For the remainder of the section, we fix c € (0, c*) and let g(s) = \pls — cs 1 / 3 for s > 0. 
The following is the key step of the proof. 

Proposition 3.2. (“No one can surf g”) For any C > 0, there exists 5 = 5(c,C ) > 0 
such that for t sufficiently large 

P {a* g N{t) s.t. Leb({s < t : X i>t (s) < g(s)}) < Ct 1 ' 3 } < e~ St1 ' 3 . 

The proof of Proposition 3.2 will take up most of this section, but we now give a brief 
justification of the result, and then show how it is used to prove Proposition 3.1. We shall 
choose a small constant /3 > 0 and let b(s) = \pls — c(s + /It) 1 / 3 for s € [0,f]. Then by 
adapting the method used in [13] for studying branching random walks, one may show 
that since c < c* if (3 is sufficiently small then for any constant K, 

P |Bi € N(t) s.t. X iit (s) > b(s) - R't 1 / 6 Vs < t} < e~ 5tl/3 

for some 5 > 0. 

Now fix K > 0 large. For large t, the function b is approximately linear on intervals 
of length Cf 1 / 3 . This will allow us to use Brownian scaling to show that if i G N(t) only 
spends time Ct 1 / 3 time below b, then it has conditional probability at least 1/2 of staying 
above b — Kt 1 / 6 , so the probability such an i G N(t) exists is also 0(e Stl/3 ). Since b < g 
this gives us Proposition 3.2. 

Before giving the details of this argument, we prove Proposition 3.1 assuming Proposi¬ 
tion 3.2. 


Proof of Proposition 3.1. We continue to write g{s) = \/2s — cs 1 / 3 , for fixed c G (0, c*) as 
above. Fix m > 0, let C = 4\/2c(l -l-m” 1 ), and let 5 = 6(c, C) be as defined in Proposition 
3.2. It suffices to show that, as t —>• oo, 

P 1 3s G [C't 1 / 3 /2,t] : d(s,m) < g(s) + 1 j —> 1. 

Next, fix t large. Recalling the notation N(t,x ) = {i G N(t) : Xi(t) > x}, let 

Ai = {* G N(t,g{t )) : Leb({s < t : X ijt (s) < g(s)}) < Ct 1/3 } 

and A 2 = {i G N(t,g(t)) : 3s G [Ct 1 ^ 3 /2,t] ■ X^s) < 0}. 


(1) 
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Also, let E be the event that d(s,m ) > g(s) + 1 for all s € [Ct 1 / 3 /2,t). On the event E, 
if i € N(t,g(t )) and i ^ A\ U A 2 then 

Mi(t) < exp(—mLeb({Ct 1 / 3 /2 < s <t : Xij(s) € (0,g(s))})) < exp(— mCt 1 ^ 3 /2). 
Since all masses are at most 1 , it follows that on E, 

y. Mi(t ) < \Ai U A 2 \ + exp(-mCt 1/3 /2)\N(t,g(t))\. 

ieN(t,g(t )) 

Also, for all y > g(t) + 1 we have C(t, 2 /) < YlieNCt g(t)) we thus have 

P jvs € [Ct 1 / 3 / 2 , t], d(s, m) > g(s) + 1 j 
= P {d(t,m) > g(t ) + 1 ,E} 


<P< y Mi(t) > m,E \ 

[ ieN(t,g(t)) ) 

< P {Ai U A 2 / 0} + P j|iV(t,g(t))| > mexp(mC't 1 // 3 /2)| . (2) 

By Proposition 3.2, for t sufficiently large, P {A± 7 ^ 0} < exp(— dt 1 ^ 3 ). Next, using a 
spinal change of measure, 

P{A 2 + 0} < E[|A 2 |] < e^ rfl/ 3 p|B(f) > - ct 1 / 3 ,3s € [Ct 1 ' 3 /2,t] : B(s) < -s/2 s} . 

Now partitioning according to the hrst interval [j, j + 1] in which B(s ) < — s/2s, 

L*J ( \ 

P jzls € [Ct 1 // 3 / 2 , t] : B(s) < —\/2s| < P < sup B(s) > s/2j > 

j=[Ct 1 / 3 /2\ U e bb + 1 ] J 

< y p( sup B(s)>s/2j\ 

j=[Ct 1 / 3 /2\ l s - J+1 J 

L*J 

= 2 y p{B(j + l)>V 2 j] 

j=\_Ct 1 / 3 /2\ 

[t\ 

<2 y exp(—j 2 /(j + 1 )) 

[Ct V3/2J 

< 2 exp(— Ct 1 / 3 /S) 


where the equality in the third line follows by the reflection principle, and the final in¬ 
equality holds for t sufficiently large. Since C > 4\/2c, it follows that 

P {A 2 / 0} < e ^ rfl/ 3 2e - ctl/3 / 3 < 2e-^ rfl/3 / 3 . 


Finally, by another spinal change of measure, 
p {\ N (t,9(t))\ > x} < x _ 1 E[|AT(t, 5 (t))|] = x _1 E 


' e -s/2 B(t) 


1 [B(t)>-ct 1 / 3 ] 


< x- l e^ ctl/ \ 


Combining the bounds on P {A\ / 0}, P {A 2 / 0}, and P {|iV(t, <?(i))| > x} with (2), we 
obtain that 


P |Vs € [Ct 1 / 3 /2,t\,d(s,m) > g(s ) + 1 j 


< e~ 5tl/3 + 2e-^ rfl/3 / 3 + m -^ e (V2c-mC/2)t^ ; 


which tends to 0 as t —>• 00 since C > 2s/2cm 1 . This establishes (1) and completes the 
proof. □ 
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For the rest of this section we work towards the proof of Proposition 3.2. We shall need 
the following lemma. Recall that we fixed c < c*. 

Lemma 3.3. There exists f3 > 0 such that for b(s) = \/2s — c(s+/3t) 1 ^ 3 and for K > 0 ,t> 
0 both sufficiently large, there exists a function A : [0 ,t] —>• [t 1 / 4 , IFt 1 / 3 ] with A(t) < /Ft 1 / 4 
and with | A'(s)| < 1 for all s € [0, t], such that for all u G [0, t] and all x € [—/Ft 1 / 6 , A (it)), 

P | b(s) — Kt x / 6 < B(s) < b(s ) + A(s) Vs < u , B(u) > b(u ) + x} 

< exp (—u — t l ^/K + y/2(A(u) — x)). (3) 

We prove Lemma 3.3 by appealing to technical lemmas from [20], which bound the 
probability that a Brownian motion stays in a narrow tube of variable width. In order 
to verify that the results of [20] apply for some function A with the above properties, 
we adapt a technique from [13]. In [13], the existence of a function analogous to A is 
constructed as the solution of a certain integral equation. We defer the details of the proof 
to Appendix A. 

From this point on, we let f3 > 0 and b(s) be as in Lemma 3.3. We assume that t is 
sufficiently large that b is increasing on [0, oo). We now show that if K is sufficiently large, 
a Brownian motion which spends at most Ct 1 / 3 time before time t below the curve b has 
a conditional probability of at least 1/2 of staying above the curve b — /Ft 1 / 6 up to time t. 

Lemma 3.4. Let (B(s)) s > o be a Brownian motion started at 0. Then given C > 0, there 
is a constant K(C ) > 0 such that for t sufficiently large, and any measurable function 
A : [0, t] (0, oo), 

P | B(s) < b(s ) + A(s) V s < t, Leb({s < t : B(s ) < 6 (a)}) < Ct 1 / 3 } 

< 2P jft(s) - /Ft 1/6 < B(s) < b(s) + A(s) Vs < t} . 

In proving Lemma 3.4, we will use the following auxiliary result. 

Lemma 3.5. Fix non-negative real numbers ( f,i > 1) For each i > 1 let 0 < 

u < tf be either a Brownian meander or a Brownian excursion of length t, t . Then writing 
T = Xa>i U, for x > 8 T 1 / 2 we have 

P < maxmaxXj(u) > xl < e _a:2 / 16T . 

( i> 1 u<ti J 

The proof of Lemma 3.5 is deferred to the appendix. 

Proof of Lemma 3.f. Write 

E = {B(s) < b(s) + A(s) Vs < t} 

Ai = { B(s ) > 6(a) - Kt 1/& V a < t} , 

A 2 = {Leb({s < t : B(s) < 6 (s)}) < Ct 1 / 3 } . 

To prove the lemma, it suffices to show that provided K = K(C ) is sufficiently large, 
P {Af | A 2 FE} < 1 / 2 , since 

P{4in/}>P{dinA 2 n/} = P{A 2 n/}(i-P{4) | i 2 n/}). 

Fix L G (Ct 1 / 3 , 2Ct x / 3 ] so that n := t/L is integer; this is possible for t large enough. 
Then, for each i < n — 2 let 6 * : [iL, {i + 2 )L\ -> K be defined by 

k(s) = b (n) + +m - b(n)) - 1 . ( 4 ) 
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Note that by convexity and since the linear terms cancel, for all i < n — 2 and s € 
[iL, (i + 2)L\, 

b(iL) H 2L (^((* 


which is less than b(s) + 1 for t sufficiently large. It follows that for t sufficiently large, 
bi < b on the interval [iL, (i + 2 )L\, for all i < n — 2. 

Next, for i < n — 2 let gi = inf{s > iL : B(s ) > bi(s)}. Also, for i < n — 2 let 
cLi = sup{s < (i + 2 )L : B(s) > bi(s)}, and let d n -2 = t. Then write 

Ui = {s € [gi, di] : B(s ) < 6 ,(s)}. 

For i < n — 2, this is the set of times when B is performing an excursion below bi which 
starts at or after time iL and ends at or before time (i + 2 )L. For i = n — 2 we have 
(i + 2 )L = t, and in this case we include a final excursion below 6 j which does not end 
before time t if it starts at or after time iL. The set U, is a union of closed intervals, which 
we enumerate as {[k.j, > 1 } according to a fixed rule (in decreasing order of size, 

say). 

For all i < n — 2, conditional on Ui, for each j > 1 the function 

(bi(li,j + s) - B(lij + s), 0 < s < r itj - l if j) (5) 

is a Brownian excursion of length r,j — lij. The case i = n — 2 is very slightly different, and 
we now describe it; for the remainder of the paragraph set i = n — 2. If B(t) > b n _ 2 {t) then 
there is no change. However, if B(t) < fr n _ 2 (f) then there there is a unique integer j > 1 
with [h j,ri j\ with = t; for this j the process described by (5) is a Brownian meander 
of length ri j — li j; for all other j the process is a Brownian excursion. All this is true even 
if we additionally condition on A 2 H E, since letting U = Ui<„_ 2 ^;, the occurrence of the 
event A 2 H E is determined by Leb {U) and By Lemma 3.5, it follows that 

P < sup( 6 j(s) — B(s)) > x 
l sEUi 

=P < sup sup (bi(s) — B(s)) > x 

[ J>1 s€[k,i,r it j] 

~ 6XP (“laLebM)) + 1 ^ <e4L ' b ( U ‘)l' (6) 

We next analyze the event A\ n A%. Note that bi + 1 is the linear interpolation of b on the 
interval [iL, (i + 2 )L\. Since b is convex it follows that b < bi + 1 on this interval. 

If A\ occurs then there is s < t such that B(s) < b(s) — Kt 1 ^. For such s, for any i 
with s € [iL, (i + 2)L], the preceding paragraph then implies that B(s ) < b(s) — Kt 1//6 < 
bi(s) - (Kt 1 / 6 - 1). 

Next suppose A 2 occurs, and suppose s < t is such that B(s) < b(s ) — Kt 1 / 6 . Then 
s is in an excursion of B(s) below b(s). Temporarily write [g,d] for the time interval 
during which this excursion takes place. Since A 2 occurs, [g, min(d, t)] has length at most 
Ct 1 / 3 so is strictly contained within in an interval [iL, (i + 2)L] for some i < n — 2. Since 
s € [g, min(d, t)] and B(g ) = b(g) > bi(g) and either d > t or B(d) = b(d) > bi(d ) but 
B(s) < bi(s), it follows that s € Ui. On the other hand, each point s lies in at most three 




- b(iL)) < b(s ) + 2L[b'(i + 2 L) - b'(iL)\ 


< b(s ) + (2Ly 


2c 


9(/3t + iL) 5 / 3 


< b(s) + 


32cC 2 
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distinct sets Ui , so on A 2 we have 

Y Leb (Ui) < 3 Ct 1 / 3 . 

i<n—2 

Finally, suppose A1nA 2 occurs. Then the observations of the preceding three paragraphs 
imply that there exists i < n — 2 and s £ Ui such that B(s) < bi(s ) — ( Kt 3 / 6 — 1) < 
bi(s) — AT 1 / 6 /2, the last inequality holding for t large. Combined with (6), this yields 

p {Ai I a 2 n E} 


<P < sup sup {bi(s) — B(s)) > Kt 1 / 6 /2, Y'' Leb(ZYj) < 3 Ct 1 / 3 

i<n—2 sEiYj 


i<n—2 


n—2 


< sup ...E' sup( 6 ,.(s) — A(s)) > Kt 1 / 6 /2 

«i+...+w„_2<3Ct 1 / 3 j = x l 


yf 2 n A 
Leb(ZYj) = Ui, A 2 n A 


Wi>0 


n—2 


< sup . E exp (—AT 2 1 1 / 3 /64tij), 

wi+...+n n -2<3Ct 1 / 3 j=l 

Wi>0 


the last bound holding provided that A " 2 > 768C so that (Kt 1 / 6 / 2 ) 2 > 64(3C't 1 / 3 ). Finally, 
letting x = A 2 t 1 / 3 /64, the function /(a) = e -x / a is convex for a £ [0, x/2], and /(0) = 0, 
so if K 2 > 384 C then for each i, f(m) < (ui/ ^ Uk) f (Y u k)■ Hence 

P {Al | A 2 CA}< e ~K 2 /192C < e -4 < 1 / 2) 


as required. 


□ 


We next state a variant of Lemma 3.4 which is proved in a similar way. 

Lemma 3.6. Let (B(s)) s >o be a Brownian motion started at 0. Then given C > 0, there 
is a constant K = K(C) such that for t sufficiently large, and any measurable function 
A : [0, t] —t (0, 00 ), u < t and z £ [b(u),b(u) + A (u)), we have 

P | B(s) < b(s ) + A(s) V s < u, Leb({s < u : B(s) < b(s)}) < Ct 1 / 3 , B(u ) > z| 

< 2P | b(s) — Kt 1 / 6 < B(s ) < b(s) + A(s) Vs < u , B(u) > . 

and 

P |a(s) < b(s) + A(s)Vs < u,Leb({s < u : B(s) < b(s)}) < Ct 1 / 3 ^ 

<2P| b(s)-Kt 1/6 <B(s ) < b(s) + A(s)Vs < u} . 

Proof. These bounds are proved in the same way as Lemma 3.4, by only considering the 
times (Ui)i when B is performing an excursion below bi on the interval [0, vi\, and using 
that z > b(u ), conditioning on B(u) > z does not affect the distribution of B on (Ui)i 
given (Ui)i. We omit the details. □ 

We are now in a position to complete the proof of Proposition 3.2, using Lemmas 3.3, 
3.4 and 3.6. 


Proof of Proposition 3.2. Choose /3 such that that Lemma 3.3 applies, and recall that 
b(s) = \/2 s — c(s + fit) 1 / 3 < g(s ) Vs > 0. Then 

P jzli € N(t) s.t. Leb({s < t : Xi jt (s) < g(s)}) < Ct 1 / 3 ^ 

< P jzli € N(t) s.t. Leb({s < t : Xj jt (s) < b(s)}) < Ct 1 / 3 ^ . 
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We shall prove that 

P {a* € N(t) s.t. Leb({s < t : X i>t {s ) < b(s)}) < Cf 1/3 } < e~ 5tl/3 (7) 

for some S > 0 for t sufficiently large, which establishes the proposition. Take K and t 
sufficiently large that Lemmas 3.3, 3.4 and 3.6 hold. Then let A : [0, t] —> [t 1 / 4 , AT 1 / 3 ] be 
as in Lemma 3.3, and in particular satisfying that A(t) < AT 1 / 4 and |A 7 (s)| < 1 for all 
s € [ 0 , t]. 

Since |A'(s)| < 1 for all s € [0, t], inf ue y J+1 ] A (u) > A (j ) — 1 for j € [0,t — 1]. Hence 
if for some i € N(j + 1 ), X t _ J+ i (s) > b(s ) + A(s) for some s € [j, j + 1], then since b is 
increasing, 

x i,j+i( s ) > Kj) + r inf . A 0 ) > & 0 ‘) + A(j) - 1 . ( 8 ) 

uepj+i] 

Using ( 8 ), and partitioning the event {3i € A(t) s.t. Leb({s < t : X^s) < b(s)}) < 
Ct 1 / 3 } according to the interval [j. j + 1] in which X t t(s) first exceeds b + A, we have that 

P jzli € N(t) s.t. Leb({s < t : Xij(s) < b(s)}) < CT 1 / 3 j 

< P 1 3* € N(t) s.t. Xi jt (s) < b(s) + A (s)Vs < f,Leb({s < t : Xi )t (s ) < 6 (s)}) < CV 1 / 3 j 

LiJ 

+E p f 3 € N{j + 1) s.t. Ajj_|_i(s) < b(s) + A(s)Vs < j, 
i=o 

Leb({s < j : A ij+ i(s) < b(s)}) < CV 1/3 , sup X ij+ i(s) > b(j) + A(j) - 1} 

ae[?jJ+l] 

< e*P ji?(s) < b(s) + A(s) V s < t, Leb({s < t : B(s) < b(s)}) < CV 1 / 3 j 

L*J 

+ ^ C +1 V{B{s) < b(s) + A(s) V s < j, 

3=0 

Leb({s < j : B(s) < b{s)}) < Ct 1/3 , sup B(s) > b{j) + A (j) - 1 }, 

where the last inequality follows by Markov’s inequality and the many-to-one lemma. By 
partitioning according to the value of A(j), we further have 


P \ B(s ) < b(s) + A(s) Vs < j, Leb({s < j : B(s) < b(s)}) < Ct 1 / 3 , sup B(s) > b(j) + A (j) — 1 
I -sebd+i] 


< P 


{ B(s) < b(s) + A(s) Vs < j, Leb({s < j : B(s ) < b(s)}) < Ct 1/3 ,B(j) > b(j) + A (j) - ^ 1/4 | 
+ P | sup B{u) > > P |b(s) < b(s) + A(s) Vs < j, Leb({s < j : B(s) < b(s)}) < Ct 1 / 3 ^ . 


Since A (j) — ^t 1 / 4 > 0, we can now apply Lemma 3.4 and Lemma 3.6 to conclude that 
P jdi € N(t) s.t. Leb({s < t : Xij(s) < b(s)}) < Ct 1 / 3 ^ 

< 2e*P { b(s) - AT 1 / 6 < B(s) < b{s) + A(s) Vs < t} 

L*J , 

+ 2^e/ +1 ( P | 6 (s) - Kt 1 ' 6 < B(s ) < b(s) + A(s) Vs < j,B(j ) > 6 (j) + A(j) - if 1 / 4 } 

3=0 V 

+ 2 e"l tl/2 p - AT 1 / 6 < A(s) < 6 (s) + A(s) Vs < j} ) • (9) 
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We can now apply Lemma 3.3 to each term. First, by Lemma 3.3 applied with u = t and 
x = —Kt 1 / 6 , since A (t) < Kt 1 ^ 4 , 

P | b(s) - Kt l/6 < B(s) < b(s) + A (a) V s < f } < exp(-i - t 1/3 /K + y/2K{t l/4 + t 1/6 )). 
By Lemma 3.3 applied with u = j and x = A(j) — ^t 1 / 4 , 

P | b{s) - Kt l/6 < B(s) < b(s ) + A(s) V s < j,B(j) > b(j) + A(j) - ±t 1/4 } 

< exp(—j - f 1/3 /AT + \/ 2 it 1/4 ). 


Finally, by Lemma 3.3 applied with u = j and x = —Kt 1 / 6 , since A(j) < Kt 1 / 3 , 

P - Kt 1 / 6 < B(s ) < b(s) + A(s) Vs < j} < exp(-j - t 1 / 3 /^ + y/2K(t 1/3 + t 1/6 )). 


Putting everything together in (9), 

P jzli G A^(t) s.t. Leb({s < t : Xi jt (s) < b(s)}) < Ct 1 ^ 3 | 
< 2 exp(-t 1/3 /iV + y/2K{t l/A + t 1/6 )) 


Ld / \ 

+ 2e ( exp(-t 1 / 3 /A" + y/2\t l/A ) + 2 exp(-t 1 / 2 /8 + 0(t 1/3 )) J 


< e 


i =o 

-5d/ 3 


for some b > 0 for t sufficiently large, which proves (7). 


□ 


4. The greatest overall particle density 

Before moving to the lower bound, we first prove logarithmic upper bounds on how the 
greatest particle density grows over time; these are needed to ensure that particle masses 
cannot decay too quickly. This may seem contradictory, but the point is that a particle 
may a priori quickly lose a large amount of mass if it finds itself in an extremely dense 
environment. The next proposition rules this out. 

Proposition 4.1. Let Z = 2 • 10 8 ; then for all s sufficiently large, 

P {sup{£(t, x) : 0 < t < s,x G M} > Zlogs} < s~ 4 . 

Proving Proposition 4.1 turns out to be a fair amount of work. In order that the idea 
is not obscured by detail, however, we set up the heart of the argument right away. 

Let z(t,x ) = Yl{i- \Xi(t)-x\<\/2} The differences between z and £ are that z only 

counts mass within distance 1/2 of x, and does not ignore the mass of particles at x (should 
there be any). 

Let z(t) = sup X z(t,x), and define a sequence (r,;, i > 0) of stopping times as follows. 
Fix s large and for the remainder of the section write N = N(s) = KLlogs. Let To = 
inf {t : z(t) > N — 1}, and for k > 0 let r^+i = inf{f > t& + 10 5 /1V : z(t) > N — 1}. Then 
i~k > 10 5 k/N, so with I = I(s) = inf {A : tj- > s}, we have I < |"As/10 5 ] and 

sup{z(t),t < s} < sup{ 2 (t),t < Tj}. 

Notice that the sequence of stopping times “ignores” small time intervals [rfe, Tfc+10 5 /A]. 
However in any time interval [r^ + 10 5 /A, r^+i), the function 2 nowhere exceeds N by the 
definition of the stopping time Tk+\- We thus have 

supj-s/i), t < s} < sup{ 2 (t),t < tj} < max A, sup sup z(t) (10) 

\ k<1 t£[i~k,Tk+W 5 /N] J 
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We prove the proposition by establishing the following facts. The first fact says that for 
k < Ns/10 5 , if z{r k ) is not too large then with high probability z(t) is not too large for any 
t € [ifc, Tfe + 10 5 /1V]. The second says that for such k, with high probability z(r k + 10 5 /1V) 
is small. 

Fact 4.2. For s sufficiently large, for all 0 < k < Ns/10 5 , 

P {sup{ 2 :(f),f € [Tfc,T fc + 10 5 /iV]} > ION, z(t/c) < N, k < /} < s~ 6 . 

Fact 4.3. For s sufficiently large, for all 0 < k < Ns/10 5 , 

P {z(r k + 10 5 /N) >N-l,k<l}< s" 6 . 

Assuming these two facts, the proposition follows easily. 

Proof of Proposition f.l. Fix k < Ns/10 5 . Note that if z(jk- i + 10 5 /1V) < N — 1 then 
z{tZ ) < N — 1. Since mass only increases by branching, it follows that almost surely a 
single branching event at time r k causes z to increase above IV—1. As all masses are at most 
1 and branching is binary, it follows that in this case almost surely z(r k ) < z(tF) + 1 < N. 
With Fact 4.3, this implies that 

P { z(r k ) > JV, k < 1} < P {z(r k _ 1 + 10 5 /iV) > N — 1, k < /} 

< P {z{r k -1 + 10 5 /N) > N -l,k-l < 1} 


We now use that for any events A, B , C we have P{dnC}<P{4nBn C}+P { B c n C}. 
By Fact 4.2 and the preceding bound, we obtain that for 0 < k < Ns/10 5 , 

P |sup{z(t), t € [ T k ,r k + 10 5 /A1]} > ION, k < /} < 2s ^ 6 

A union bound and (10) then yield 


t<S 


P < sup z(t) > 10A^ > < P < sup sup z(t) > ION 


k<I te[r i .,r fc +10 5 /N] 
LATs/10 5 J 

£ E p 


k =0 


sup z{t) > ION, k < I 

t&[T k ,T k + 10 5 /N] 


Sll+ ^)' 2s " 


< s 


-4 


the last inequality holding for s large. Finally, it is easy to see that sup x ((t, x) < 2 z(t), so 
the same bound holds for P {sup t<s sup T f(t, x ) > 201V}, which proves the proposition. □ 


The reader who is willing to believe the Facts 4.2 and 4.3 without proof - or who is 
impatient to see how Proposition 4.1 is used to prove the lower bound from the main 
theorem - could skip directly to Section 5 at this point. 


4.1. Proofs of Facts 4.2 and 4.3. We first prove a handful of technical estimates re¬ 
quired for the proofs. The first shows that a fixed mass of particles is extremely unlikely 
to quickly increase its total mass. Recall the definition of P x ,m from just before the start 
of Section 3. 

Lemma 4.4. Fix x = (x \,..., x k ) € and m = (mi,... , m k ) € [0,1] A: . Under P x , m , 
for 1 < j < k let Gj(s) = ff{i : jj, iS (0) = j} be the number of time-s descendants of Xj. 
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Then for any Jc {I,..., k}, any x > Yhjej m j> f or aR ^ — log 2 and all 5 > 0, 


YmjGjit) > (1 + 5)x > < 2 ( 2 ^(l-e-y) : 


, j'eJ 


Proof. We may clearly assume J = {1,... ,k}. Also, adding particles to increase the mass 
of the starting configuration can only increase the probability we aim to bound, so we 
may assume that x = The random variables (Gj(s),l < j < k) are i.i.d. and 

are Geom(e _s )-distributed (see, e.g., [18]). Lemma B.2 provides upper tail bounds for 
weighted sums of geometric random variables where the individual coefficients are small 
compared with their sum. Using that lemma (with e = 1 — - this is where we require 

that t < log 2 ), the result follows. □ 


Since Gj(s) is non-decreasing in s, we have 

sup Mi )S {0 ) = sup rrijGj(s) = 'Y m jGj(t ). 

se[M Uii,A 0 )eJ} «e[o,t] i6J j&J 

Combining this with the preceding lemma thus also yields the following bound. 

Corollary 4.5. With the hypotheses and notation of Lemma 4-4> 

P x ,m < sup Y M ijS (0) > (1 + S)x 

l se[0 ’ t] B:L,a(0)eJ} 

The next proposition says that mass does not travel far in a short time, even once 
branching is taken into account. 

Proposition 4.6. Fix x = (aq,..., aq) € and m = (mi,..., m*,) € (0, l] fc . Then for 
all x > Yli<i<k m i> f or dlt>0,L> 1/2 and v > 0, we have 

„ . j exp(i — L 2 /(2t)) 

Y M i)t ( 0 ) > vx ^ ^- /V ■ 

{v.Xi(t)-X itt (0)>L} ) 

Proof. We begin with a few simplifying assumptions. First, we may clearly assume that 
Xi = 0 for all i < k. Next, adding particles to the system at time 0 can only increase the 
probability we aim to bound, so we may assume that x = Yli=i m i- 

For j < k write Sj = {i < n{t ) : ji t(0) = j} for the set of indices of time-t descendants 
of Xj. Then let Rj = {i € Sj : Xi(t ) — A,;^(0) > L}, so that 

k k 

Y M i, t ( o) = YY1 m m(o) = Y m j\ R j \ ■ 

{i:Xi(t)-X itt (p)>L} j=li£Rj j=l 

By the many-to-one lemma, for W a one-dimensional Brownian motion, 




E [|i^|] = e*P {W t -W 0 > L}< exp(f - L 2 /(2t)) 
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for L > 1/2. This bound does not depend on j < k. It then follows by Markov’s inequality 
that for v > 0, 

Mij (0) > vx > = P < nrij\Rj\ > vx 


{i:Xi(t)-Xi, t (0)>L} 


E 


< 

< 


j<k 

Ej<k m j\ R j\ 


VX 

exp(f — L 2 / (2t)) 
v 


□ 


where we have used in the last inequality that J2j<k m i = x - 

In the sequel, we in fact use the following corollary, which extends Proposition 4.6 by 
considering all times in an interval [0, t], rather than a fixed time t > 0, at the cost of a 
slightly weaker bound. 

Corollary 4.7. Under the conditions of Proposition f.6, for all to > 0, L > 1/2 and 
v > 0, and all x > Yli<k 


x,m ^ sup yy M ijt ( 0) > 2vx } < 

{i:Xi(t)-X i}t (0)>L} 

Proof. Consider the stopping time 


2 exp(t 0 - L 2 / (2f Q )) 


r = inf < t 


E 


Mij( 0) > 2vx 


{i:Xi(t)-X i:t (0)>L} 


By symmetry, 


yy Mi tto ( 0)>vx 

{i:Xi(to)—Xi : t 0 (0)>L} 

and the corollary follows. 


r < t 0 }> 


2 ’ 


□ 


The next lemma says that a large, concentrated mass will quickly decay; once we prove 
this we will have all the tools we need to establish Facts 4.2 and 4.3. 

Lemma 4.8. There exist to > 0 and C > 0 such that the following holds. Fix x = 
(xi, ... ,Xk) € and m = (mi,..., m*,) € [0, l] fc . Let J = {j : \xj\ < 1/4}, and suppose 
A = > C. Then for all t € [500/4., to], setting I = {i : ji,t(0) € J} we have 


yy M t (t) > 4/24 i < 2e~ 200A 
. iei J 


Proof. The proof is divided as follows. First, the total mass at time t of particles whose 
trajectory branches at least once is small. Next, among non-branching trajectories, the 
total mass which moves far from the origin is small. Finally, particles whose trajectories 
do not branch and stay near the origin will lose a large amount of mass since they are a 
dense environment. We now formalize this. 

Write = {* G I : / i,ji,t( 0) = jV(0)} for the indices of particles starting 

near (distance < 1/4) to the origin whose trajectories branch before time t. Then let 
I \ /ft = If U I n , where 

If = \iel\l b : |X M (0)| < 1/4, sup |X M (s)[ > 1/2 i 
I s6[0,t] I 
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indexes non-branching trajectories that start near the origin but move far (distance >1/2) 
from the origin before time t, and where I n = I\(IfUl b ) indexes non-branching trajectories 
that stay near the origin. Then with M b = Y2iei b Mi(t) and Mf, M n defined accordingly, 
we have 

Y Mi(t) = M b + M f + M n . 
iei 

We begin by considering branching trajectories. For each 1 < j < k, let Gj = #{i € 
I '■ ji,t( 0) = j}■ Then i € /& precisely if (0) € J and G Jit (o) > 1- Since masses decrease 
with time, 

Y Mi ® <E^, t (0) 

i£lb 

= J2 m i G i 1 [G j > i]- 

3&J 

Next, since the Gj are integer-valued, 

^nrijGjl [G .>i] = Y m j( G j - 1) + J] mjl[ Gj >i] < 2 ^ m A G i ~ !)» 
jeJ j&J j&J j£J 

which with the preceding bound gives 

Mi(t ) < mjGj — A). 

iei b jeJ 

By Lemma 4.4, it follows that for any fixed 5 >0, if t < log 2, 

P I M ^ f ) > 26A \ < p < Y m i G i - ( x + 


ie4 


J£ J 

< 2(2 1+<5 (1 - e" t )' 5 ) A 

< 2(2 1+< V) A 


<e~ 200A , (11) 

the last bound holding for t sufficiently small that 2 2+s t s < e~ 200 . We next bound 
Z)iS/n the total final mass from “typical” trajectories, which do not branch and 

do not move far from their starting position by time t. Fix c € (0, 1) and let E be the 
event that for all s G [0, t], Yl{i-jXi{s)\<i/ 2 } Mi(s) > cA. On E. if i g I n has jij( 0) = j then 
Mj(t) < rrij ■ e~ tcA . We thus have 

Y Mi(t) 1[ B] <Y m r e ~ tcA ' ![E] = Ae~ tcA 1 [B] . 

i£ln j&J 

Next, let 7 n (s) = {j/,t(s) : i € 7 n } be the indices of time-s ancestors of individuals in 7 n . 
Since trajectories indexed by I n do not branch, X)ie/„(s) Mi(s) is decreasing for s € [0,7]. 
Necessarily |7Q(s)| < 1/2 for i € 7 n (s), so if E c occurs then there is s € [0,7] such that 
Y.ieUs) Mi( s ) — cA We thus have 

yy ^7,;(1)1 [ec] < cy4l[ E c], 

and the two preceding bounds together give 

yy Mj{t,) < max (cA, Ae~ tcA ^j . (12) 

Finally, we turn to the final mass of non-branching trajectories that move far from the 
origin, counted by Yliei f For any i G 7, If j/,t(0) = j and \xj\ < 1/4 then in order to 
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have sup se r 0 ^i |Xj ;t (s)| >1/2 the trajectory leading to Xi(t) wanders a distance of at least 
1/4 from its starting position. Let W denote one-dimensional Brownian motion started 
from the origin. By the reflection principle and the fact that P {G > x} < e~ x / 2 /(2a:) for 
G a standard normal and for all x > 0, we have 

P |sup \W S \ > 1/41 < 4p| W t > < 8exp(—1/(32£)). 

Since an individual trajectory of X has the law of Brownian motion, for a particle starting 
at distance less than 1/4 from the origin whose trajectory never branched, the above is a 
bound on the probability the trajectory attained distance 1/2 from the origin. It follows 
that 

y ] A4j(t) > s t y ] Trij ■ , 

ielf j£J 

where the terms are i.i.d. Ber(8exp(—l/(32f))). The variance of the latter sum is 
bounded by A ■ 8exp(—1 /(32t)), so Theorem B.l yields that for any fixed b > 0, 

bA 

I / O-i— LUSZt) \ 

> (b + 8exp(—l/(32f)))A 

icif 

the final inequality for t sufficiently small. 

We now combine (11), (12) and (13). This yields that for t sufficiently small, and in 
particular provided that 2 2 + s t s < e -200 , ((8e 1 ~ 1 ^ 32t ' ) )/b) b < e -200 and that 

25 + max(c, e~ tcA ) + b + 8exp(—l/(32t)) < 


< 


' 8e l—1/(32 1 )' 


< e 


-200A 


(13) 


24 


we have 


P Mi - A / 24 | < 2e~ 200A 


It can be checked that taking 5 = b = c = 1/100 does the job when t > 100 log 100/A 
(so that max(c, e~ tcA ) = 1/100) and t is sufficiently small (it is in order to satisfy these 
simultaneously that we require a lower bound on A). This completes the proof. □ 

Proof of Fact f.2. Let Z/2 = {y/2:y€Z}. Define the event 

E = {max{|Xj(r)|, i >l,0<r<s + 10 5 /IV} < 3s} . 

Any unit interval [x — 1/2, x + 1/2] is covered by at most two intervals from {[y — 1/2,y + 
1/2] : y € Z/2}. It follows that on E, if < s but sup {z(t),t € [rfc,Tfc + 10 5 /A"]} > 101V 
then there is y with y € [—3s, 3s] n Z/2 such that 

sup ^2 Mj(f) > 5 N. 

te[T k ,T k +10 5 /N] 

When k < I we have tj~ < s, so 

P < sup z(t) > ION, z{jk) < N, k < I > 

I te[T k ,T k +W 5 /N] I 


sup z(t, y) > 5IV, z(jk) < N 
te[r fe ,r A .+10 5 /A''] 


(14) 


<P{£ C }+ P 

ye[-3s,3s]n7./2 

Our bound on the above summands works identically for each y € [—3s, 3s] 0 Z/2; we 
explain it for y = 0 to avoid notational overload. So we wish to bound 


sup z(t, 0) > 5N,z(rk) < N 

t£[T k ,T k + 10 5 /N] 
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Our strategy is as follows: we use Corollary 4.5 to show that with high probability, 
for all t € [r k ,r k + 10 ,5 /1V] the total contribution to z(t, 0) from descendants of particles 
with \Xi{jk)\ < 3/2 is at most AN. We then use Proposition 4.6 to show that with high 
probability the contribution to z(t, 0) from descendants of further-off particles decreases 
quadratically (as a function of |W(t/c)|); since the quadratic series converges, this implies 
a bound on the total contribution from far-off particles. We now proceed to details. 

For n € Z let 


Y n = sup ^2 Mi(t ); 

te[r k ,r k +10 5 /N) {i .\ Xi {t)\<\/2,\X^ t {T k )-n\<l/2} 


Y n counts the greatest contribution at any time t € [r k , r k + 10 5 /IV], to the mass near 0 
from particles that at time r k are near n. We clearly have 


sup 

te[r fc ,r fc +10 5 /JV] 


Z{t, 0 )<£ Yn ■ 
nEZ 


(15) 


As sketched above, we bound the sum in two parts: the contribution from Y_i, Yq and Y\ 
is handled separately from the rest, and we do this first. Note that since masses decrease 
with time, 


Y-! +Y 0 + Yl < sup ^2 MiA T k )• 

t£[T k ,r k +10 5 /N] {i .\ Xiit (r h )\<Z/2} 


If z{jk) < N then X){r|x i (r fc )|<3/2} ^i{ T k) < 3Ai so, by Corollary 4.5 and the strong Markov 
property, 


P {F'_i + Y 0 + Yi> AN, z(r k ) < N} < 2 ^2 1+1/3 (1 - e - lo5 /A7)i/ 3 ^ 3iV 

< (20 5 /N) n . 

Now consider n € Z with |n| > 2, and assume by symmetry that n > 0. If |Xj(t)| <1/2 
but | Xit(rk) + n\ < 1/2 then Xi(t) — Xij{r k ) > n — 1. Assuming z(tx-) < IV, in particular 
we have z(r k , —n) < IV. Furthermore, 

Y- n < sup Mjt(rfe). 

t6[Tfc , rfc+105/ iv ]{ . :Xi(t) _4r (Tfc)>n _ 1} 


When n > 2, applying Corollary 4.7 with to = 10 5 /1V, L = n — 1, v = 1/(20(n — l) 2 ) and 
x = N, we then obtain that 


P P-n > 


at 1 /in 5 

, 2(r fc ) < N \ < 40(n - l) 2 exp f — 


10(n - 1) 2: 


< exp — 


N(n-lf 


3- 10 5 


C s -10(n-l) 2 


N(n — 1) 2 A 

2 • 10 5 J 


(16) 


The final inequality holds since N = N(s) = 10‘ logs; the second inequality holds provided 
N is sufficiently large. We emphasize that once N is large enough the inequality holds for 
all n > 2. Note that by symmmetry the same bound also holds for Y n . 
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Using (15) and the two preceding probability bounds (and the fact that (1/10) 'f2\ n \>2( n ~ 
l ) -2 = 7t 2 /30 < 1), we thus have 



where the last inequality holds for s sufficiently large. The same argument yields the same 
bound with z(t,y ) in place of z(t, 0), and (14) then gives 


P < sup z(t) > ION, z{jk) < N,k < I > 

[te[r fc ,r i .+10 5 /7V] J 

<P{E c } + 12s. ^00 + 4s -10 ^ 

<p{u c } + s - 8 , 

the latter bound holding for s large, since N = 10 7 logs. To conclude, we use the fact that 

P{max{|X;(s + 10 5 /N)\,i > 1} > 3s | E c } > i, 

which follows by considering the stopping time r = infjr : max{|Xj(r )|,i > 1} > 3s} and 
using symmetry. This yields 

P {E c } < 2P{max{|Xi(s + 10 5 /iV)|,i > 1} > 3s} 

< 4E [#{i : Xi(s + 10 5 /N) > 3s}] 

= 4e'P {iV(0, s + 10 5 /iV) > 3s} 

< e- 3s (18) 


where the last two inequalities hold for s sufficiently large. 


□ 


Proof of Fact f.3. The proof has aspects which will be familiar from the previous proof; 
we describe these first. We recycle the event E from the preceding proof. Note that on 
E n {k < 1} we have 

z{r k + 10 5 /N) < 2 sup z{r k + 10 5 /iV, y ), 

2 /S[—3s,3s]nZ/2 
so 


P { z(r k + 10 5 /AO > N - 1, z(Tk) < N,k < I, E) 

< ^2 P iz(r k + 10 5 /N,y) > N ^ 1 ,z(r k ) < n\ . (19) 

t/e[-3s,3s]nz/2 ^ ' 
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We once again focus on the case y = 0 for notational simplicity. We write 

= 22 Mi(r k + 10 5 /N ). 

{i:\Xi(r k +10 5 /N)\<l/2,\X iTk+lo5 /N (r k )-n\<l/2} 

The indices of summation correspond to particles with position near 0 at time r k + 10 5 /N, 
whose time r k ancestor had position near n. We have 

z{r k + 10 5 /N, 0)<£ Z n- 

nEZ 

Now similarly to the argument leading to (16), apply Corollary 4.7 with to = 10 5 /1V, 
L = n — 1, v = l/(40(n — l) 2 ) and x = N to bound Z n for \n\ > 2. We obtain that for s 
sufficiently large (since (1/20) X^|n|>2( n — 1)~ 2 = x 2 /00 < 1/2) 

P |z(r fc + 10 5 /1V,0) > N 2 1 , z(r k ) < IV| 

<P [z-i+Z Q +Zi>—,z(T k ) < nX + 22 p l z n> _ 1 N 2 ^( r fc) < N \ 

) {neZ:\n\>2} ^ 1 } 

<P |z_! + Z 0 + Z l > J,z(r k ) < nJ + 4s -10 . (20) 

We now bound Z__i + Zq + Z\ from above by the total mass at time r k + 10 5 /N of 
individuals whose time-r k ancestor lies in [—3/2, 3/2], More precisely, recall that Xj^(s) 
is the (location of) the time-s ancestor of Xj(f), and write 


D e = 22 M ^ Tk + 105 /^)- 

{^i, Tk + W 5,N( T k)e[mM+^)/m 


Then 


z~ i + z 0 + Zi < 22 Dg ■ 


te[-3,2]nz 

This holds because the time-Tfc ancestors of particles counted by Z_i + Zq + Z\ all lie 
in [—3/2, 3/2] = U^e[-3 2]nz[^/2> + l)/2]. The bound may be strict because particles 

counted by Z_i + Zq + Z\ are additionally required to lie near 0 at time r k + 10 5 /N. 
Bounding each of the summands Dp by the largest summand, we then have 


Z_ i + Zq + Z\ < 6 max Di , 
te[-3,2]nz 


so 


P |z_! + Z 0 + Zi > J, zfa) < jv| 

<6 max P { > —, z(r k ) < N \ 

.fe[-3,2]nz (24 J 


The final probabilities are not hard to bound: if Dg hearkens from a total time-Tfc mass 
which is very small then at time r k + 10 5 /N it is still rather small by Corollary 4.5. On 
the other hand, if the aggregate mass of its time-Tfc ancestors was larger (but still at most 
N) then by Lemma 4.8, at time r k + 10 5 /lV that ancestral population has lost most of its 
mass. 

More precisely, since Mi(r k + 10 5 /1V) < M i Tk+10 5/jvCnt) f° r each i, by Corollary 4.5 and 
the strong Markov property, 


P 


Di > 


N 

24 


E 


Mj(n) < N/is 


{j:X j (r k )e[i/ 2,(<+l)/2]} 


< 2(4(l-e" 1 ° 5/iV )) iV/48 < 2 


^2o 5 y v/48 





20 


LOUIGI ADDARIO-BERRY AND SARAH PENINGTON 


Now assume that N = N(s) = 10 7 logs > 48 C, where C is the constant from Lemma 4.8. 
By that lemma, since 10 5 /iV > 500/(lV/48), 


P {D t > 


N 

24 


Y Mjijk) € [N/A8,N] \ < 2e~ 200N/48 < e~ AN 

{r-Xj{T k )£{t/ 2,(<+l)/2]} 


the latter inequality for N = N(s) sufficiently large. This bound holds for each t € 
[—3,2] n Z. Under the assumption that z(r k ) < N, one of the conditions in the above 
conditional probabilities must occur. It follows that 


6 max P 
£e[-3,2]nz 


N 

D i > ^ z ( T k) < N 


< 6 max 


( ( 20 5 \ N/i8 
{ 2 {ir) 



so for N sufficiently large 

P jz_i + Z 0 + Z 1 > j,z{r k ) < jvj < 6e~ 4N = 6s" 4 ' 107 . 


Combined with (20) this gives 

P \z{r k + 10 5 /iV,0) > N 2 1 ,z{r k ) < iv| < 5s' i0 . 

The same bound holds for each z(rk+10 5 /N, y), so using (19) and the bound P {E c } < e^ 3s 
from (18), for s large we obtain 

P {z(T k + 10 5 /JV) > N - 1, z(r k ) < N,k < 1} < 60s" 9 + e" 3s < s" 8 . 

The proof is almost complete; to finish it off we need to deal with the event { z(r k ) < 
A^} in the preceding probability. To do so we use induction. First, for s large, since 
N = N(s) = 10 7 logs and tq = inf{t : z(t) > N — 1}, then z(tq) < N — 1. It follows that 
z(to) < z(t 0 ) +1 < N (this was explained in the proof of Proposition 4.1), so when k = 0 
we have 


P {z(r k + 10 5 /AO >N-l,k<l} = P {z(r k + 10 5 /N) > N - 1, z(r k ) < N,k < /} < s -8 . 

For larger k, similarly if z(r k -\ + 10 5 /A^) < N — 1 then z{r k ) < z(t j7) + 1 < N. We thus 
have 


P{z(T k + 10 5 /N) >N-l,k<l} < P {z(T k + l0 5 /N) > N — 1, z(r k ) <N,k<l} 

+ P {z(r k ) > N, k < 1} 

< s~ 8 + P {z(r k - 1 + 10 5 /iV) > N - 1, k - 1 < /} , 
so by induction and the hypothesis that k < Ns/10 5 , 

P {z(r k + 10 5 /N) > N - 1, k < 1} < {k + 1) • s" 8 < Ns~ 7 < s~ 6 . □ 


5. Lower bound 

The next proposition restates the second inequality of Theorem 1.1. Recall that c* = 

34 / 3 ^/ 3 / 27/6 _ 

Proposition 5.1. For any m € (0,1), almost surely 

lim inf 

t->-oo 


\/2 1 - D(t,m) * 

- m -< c 

t 1 / 3 
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Given a function / : [0, oo) -A R, for t > 0 let I(t, /) = {* > 1 : Vs € [0, t], Xij(s) > 
f(s)} be the indices of particles whose ancestral trajectory stays above / up to time t. 
Note that \I(t,f)\ is decreasing in t: if a trajectory stays above / to time t then it also 
stays above / to time t' < t. It follows that P {Vt, I(f, f) / 0} = Hindoo P {/(/, f) / 0}, 
and this is a decreasing limit. We will use the following result of Roberts [20]. 

Lemma 5.2 ([20], Theorem 1). Let g(t ) = \f21 — c*t 1//3 + c*t 7 ^ 3 / log 2 (t + e) — 1. Then 

lim P {I(t , g) / 0} = p* >0 

£—>■00 

The idea of the proof of Proposition 5.1 is that if the density is always low beyond g 
then a particle staying beyond g will have reasonably large mass at time t; the lemma 
guarantees that such a particle has a reasonable chance p* of existing. The next corollary 
implies that at the cost of a constant shift of the function g, we may increase p* as close 
to one as we like. For c € R write g — c for the function with ( g — c)(x) = g{x) — c. 

Corollary 5.3. Let C* = inf{c : Vt, I(t,g — c) / 0}. Then almost surely C* < oo. 


Proof. The proof technique is sometimes called an amplification argument. Consider the 
n(t) ~ e l independent copies of the BBM rooted at time-f particles, the Pth copy having 
initial individual at position Xi{t). Suppose the “translate by Xi(t) v of the event from 
Lemma 5.2 occurs in the /Pth copy; more precisely, suppose that for all t! > t there is a 
descendant Xj(t') of X k (t) such that for all s € [t, t'], 

Xj, t >(s) - X k (t) > g(s -t)> g(s) - V%t - c*t 1/3 . 

For s < t we also have 


Xj t r(s ) > inf Xi(s) > inf inf Xj(s) > g(s) + inf inf Xjfs) — supg(s), 

’ *>1 s<ti>l s<ti> 1 

so in this case 

C* < — inf inf Xi(s) + V2 1 + c*t 1//3 . 
se[o,t] i>i 

By the branching property (i.e. the independence of the trajectories emanating from each 
of the particles (Xi(t),i > 1)), it follows that 


P f C* < 3t + V2t + 


< P{n(t) < 2*) + p| inf inf XRs) < -3ti + (1 - p* 
lse[o,t]i>i I 


\2 t 


( 21 ) 

where p* is the constant from Lemma 5.2. Since n(t) is Geom(e~ 4 ) we have P {n(t) < 2*} < 
(2/e)*. Finally, let a = inf{s : inf^i X l (s) < —3t}, so inf se r 0j ji infj>i Xj(s) < —3 1 if and 
only if a < t. Considering the descendants of the first individual to reach position —3 1, by 
symmetry we have 

p{ inf X t (t) < -31 


i> 1 


cr <t> > - , 


so 


P {a < t} < 2P |inf Xiit) < -3f| < 2e*P{iV(0,t) < -3 1} < e~ 7t/2 . 
These bounds and (21) then yield 

P|C* < 3t + V2t + c*t 1/3 } < (2/e)* + e~ 7t/2 + (1 - p*) 2 * . 

This can be made arbitrarily small by taking t large. 


□ 


In order to prove Proposition 5.1, we require one final lemma which shows that a small 
mass will quickly increase to form some region of constant density within a constant 
distance. 
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Lemma 5.4. For all e > 0 and m £ (0,1) there is C > 0 such that for all x £ and 
m £ (0, l] k , if z := £{i:|n|<i} m > 0 then 

Px,m{3t € [ 0 ,C (1 + log(l/z))],x £ [-(7,(7] : C (t,x) >m}> 1 - e. ( 22 ) 

To prove the lemma we use the following fact, whose proof is left to the reader. 

Fact 5.5. For all e > 0, there are to = to(e) and c = c(e) > 0 such that 

P {Vt > t 0 , #{z : Vs £ [ 0 ,t], |X M (a)| < c} > (e - e)*} > 1 - e. (23) 

One straightforward way to prove the fact is as follows. First show that p(c), the survival 
probability of branching Brownian motion with absorbing boundaries at —c and c, started 
from the origin, satisfies p(c) —>• 1 as c ^ oo. Then use a suitable branching process 
approximation. As an aside, we note the very nice recent work [ 11 ] on the asymptotics of 
this survival probability for c near the critical width c below which p(c) = 0 . 

Proof of Lemma 5.f . The claim is clearly true if z > m, and we hereafter assume z £ 
(0, m). We also assume e is small enough that (e — e)e _m (l — e 1 / 2 ) > (1 + e); this can only 
make our job harder. 

By relabelling, we may assume that for some 1 < k! < k we have \xf < 1 for 1 < z < k! 

and |xj| > 1 for z > k', so that z = Yli<i<k' m *- We a ^ so assume x h ■ • •, Xk are ordered so 

that (mi, 1 < z < k') is decreasing. 

For 1 < z < k' let Ji(t) index the time-t descendants of x' ? ; whose trajectory stays fairly 
near the origin, i.e., 

Ji(t) = > 1 : je,t( 0) = i, \X ttt ( 3 ) - Xi\ < c Vs £ [ 0 , t}} , 

where c is chosen as in Fact 5.5. By that fact, we then have for to = ^o( e ) 

P {Vt > t 0 , | J»(t)| > (e - e)*} > 1 - e. 

For 1 < n < k' write 

S n = #{1 < i < n : Vt > t 0 (e), |«7*(t)| > (e - e) 4 }. 

Then S = (S n )i< n <k' stochastically dominates a random walk with Bernoulli(l — e) steps. 
It follows by a ballot-type theorem ([14, Corollary 11.17], for example, is sufficient for our 
needs) that for any A > 1, 

P{3rz<fc' : S n < (1 — Ae)rz} < A" 1 . (24) 

We hereafter assume t > to(e). Now suppose that ((s, x) < m for all s < t and |x| < c+1. 
Then for each 1 < i < k', for all j £ J*(t), Mj(t) > mt ■ e~ mt , so 

X) X M j(t)>e- mt rn i \J i (t)\>e- mt (e-e) t Y m * 1 [NiWI>(e-e)‘] • 

1 <i<k' jeJi(t) 1 <i<k’ 1 <i<k' 

Since the masses m; are decreasing in z £ { 1 , 2 ,... , k'}, if S n > (1 — e 1 / 2 )n for all n < k' 
then it follows that 

l<i<k'jeJi(t) i =1 

By our assumption on e, we have e~ m (e — e) > (1 + e)/(l — e 1 / 2 ) > 1 + 2e, so this gives 

Y Mj(t) > (1 + 2e) t ~ 1 z > c + 2 , 
j:\Xj(t)\<c+l 

the last inequality provided that t > 1 + log 1+2e ((c + 2 )/z). Since [—c — 1, c + 1] can be 
covered by [c + 2 ] intervals of radius 1 , we see that in this case there is x with |x| < c + 1 
such that C,(t,x) > 1 . 
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To sum up: assuming the random walk S behaves, and that t > to(e) and t > 1 + 
l°Si+ 2 e(( c + 2)/z), either £(s,®) > m for some s < t and \x\ < c+ 1, or else £ (t,x) > 1 for 
some x with \x\ < c+1. By taking A = e -1 / 2 in (24) and choosing C = (7(e) appropriately, 
we obtain 

Px,m{3s € [0,C(1 + log (1/+))],£ € [-<7,(7] : ({t,x) >m}> 1 - e 1/2 . □ 

We are now ready for the final proof of the paper. 

Proof of Proposition 5.1. Fix m € (0,1). Let Z = 2 • 10 8 and 

t* = inf{r > 0 : Vt > r, sup sup£(s, x) < Zlogt}, 

se[o,t] xeR 

and note that t* < oo almost surely by Proposition 4.1 and the first Borel-Cantelli lemma. 

Fix e > 0 and by Corollary 5.3 choose L > 1 large enough that P {max(C*, t*) > L} < e. 
Fix t much larger than L (so that log log t > L, say). 

Let a = inf{s > t 1 / 4 : D(s, 1 /t) > g(s) — (7* — 1}. We first suppose that a > t, so that for 
all s € [t 1 / 4 , i] we have D(s , 1/t) < g(s) — C* — 1. Let i* be such that Xi*j(s) > g(s) — C* 
for all s € [0, t}] such i* exists by the definition of C*. If t* < L < t then we have 

— log Mi* (t) = [ ((s,Xi* it (s))ds 
Jo 

P 1 ' 4 P 1 

< / ((s,Xi*j(s))ds + / -ds 

Jo Jt 1 / 4 4 

< Zt 1//4 log t + 1 , 

the last bound because when t > t* the integrand is at most Zlogf. 

Let (7 = C(e,m) be the constant from Lemma 5.4. Then by that lemma (applied with 
z = Mi* (t) > exp(—1 — Zt 1 / 4 logt)) and the Markov property, given that {t* < t}, with 
probability at least 1 — e there is s € (t,t + (7(2 + Zt l/4 logt)) and x with |x| < C such 
that £(s,Xj*(t) + x) > m. If this occurs, and additionally C* < L we have 

D(s, m) > Xi* (t) — C > g(t ) — C* — C > g(s) — s 1 / 4 log 2 s , 

the last bound holding for all t sufficiently large since s — t < (7(2 + Zt 1 / 4 log t), and for s 
and t large we have g(s) — g(t ) = 0(s — t). We thus have 

P jzls > t : D(s, m ) > g(s) — s 1//4 log 2 s \ a > t j 

>P |max((7*, t*) < L, 3s > t : D(s, m) > g(s) — s 1 / 4 log 2 s \ a > t j 

— P {max((7*, t*) > L \ a > t} 

>1 — e — P {max((7*, t*) > L \ a > t} . (25) 

Next suppose that a < t. Apply the strong Markov property at time a, and apply 
Lemma 5.4 just as above (but with a starting mass in [D(a, 1/t ) — 1 ,D(a, 1/t ) + 1] of at 
least 1/t = e~ logt rather than 4lo § t ). We obtain that with probability at least 

1 — e there is s € (a, a + (7(1 + logt)) such that 

D(s, m ) > g(cr) - C - C* > g(s) - log 2 s , 

the last bound holding for t sufficiently large since s — a < (7(1 +logt) and logt < 4 log a < 
4logs, and under the assumption C* < L. 

Since a > t 1 / 4 and log 2 s < s 1 / 4 log 2 s, it follows that 

P |^s > t 1 / 4 : D(s, m) > g(s) — s 1 / 4 log 2 s \ a < 1 1 

>1 — e — P {max((7*, t*) > L \ a < t} . 
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Now combine this with (25) using the law of total probability. We chose L large enough 
that P {ma x(C*,t*) > L} < e, so we obtain 

P 13s > t 1 / 4 : D(s, m ) > g(s) — s 1 / 4 log 2 s| > 1 — 2 e . 

Finally, if D(s,m) > g(s) — s 1 // 4 log 2 s then 

y/2s — D(s, m) * c* 1 log 2 s 

s 1 / 3 — log 2 (s + e) s 1 / 3 s 1 / 12 ’ 

which tends to c* as s —>• oo. □ 


6. Discussion and questions 

• The analysis of the paper should carry through fairly straightforwardly to higher 
dimensions R fc , provided we redefine d(t,m) and D(t,m) as 

d(t,m ) = min(|.x'| : £(£,x) < m}, D(t,m ) = max{|.x'| : ((t,x) > m} . 

At time t, the density is then at least m within the ball of radius d(t, m) around 
0, and less than m outside the ball of radius D(t,m ) around 0. The proof of the 
lower bound is then the same as in Sections 4 and 5. The proof of the upper bound 
requires ruling out the possibility that the modulus of a particle in the BBM stays 
ahead of a moving barrier g even though it cannot have consistent displacement 
more than g in any fixed direction. In order for our proof techniques to carry 
over, this requires sample path estimates for Bes(/c) processes analogous to the 
ones derived in this work for Brownian motion. We expect such estimates to hold 
for all k > 1 , though verifying this may be technical. 

• We believe that Proposition 5.1 predicts the “true” front location, in that both 
D(t,m ) and d(t,m) are typically at distance o(t 4 / 3 ) from \/2 1 — c*t 4 ^ 3 when t is 
large. This is our justification for the remark in the final paragraph of Section 1. 

• In the same way as the KPP equation describes the evolution of multiplicative 
functionals of BBM [18], it seems plausible that the model proposed in this work 
(or a related model) should be connected to an equation of the form 

u t = \u xx -u(l-u)~ I u(t, y)dy. 

z d {y-\y-x\<i} 

This equation has steady states at 0 (unstable) and 1/2 (stable), and is redolent 
of a family of “non-local” Fisher-KPP-type equations which was introduced [5] to 
model populations in which aggregation can have both a competitive advantage 
(safety in numbers) and disadvantage (due to competition for resources). These 
equations have received substantial study [2, 6 , 9]; the survey [22] contains many 
further references, as well as perspective on the biological motivations for such 
study. 

If a probabilistic model for such an equation were found, it could yield new re¬ 
sults on, e.g., the front propagation speed or temporal fluctuations of solutions to 
the above equation. Conversely, a glance at that literature suggests new probabilis¬ 
tic questions: for example, what if the effect of competition is described by a kernel 
k, where k(\x — y |) describes the degree of competition for resources between indi¬ 
viduals at spatial positions x and yl In our model we took n(\x—y\) = l[[ 3 ; _ y | e (o i i)]; 
a kernel which allows substantial long-range interaction might yield rather different 
dynamics. 

• As mentioned in the introduction, one may reasonably consider the mechanism 
for mass growth in our model - both children inherit the mass of the parent - 
nonphysical. More physically realistic (at least for amoebae) is for the children to 
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each have half the mass of the parent. One must also then change the rules to 
allow for mass growth; a reasonable modification is to take 

C (t,x)= 

and 

Mi(t) = exp( f (1 - ((s,Xi(s)))ds ). 

Jo 

In other words, the mass of an individual can increase, when there is little nearby 
competition for resources - but the larger particles get, the harder it is for them 
to sustain themselves. The key point is that 1 is still a universal upper bound on 
the greatest mass of any particle. 

We conjecture that any lack of physical realism in our model is relatively in¬ 
significant for the long term behaviour, and more concretely that the front location 
behaves similarly in the two models. As partial evidence for this, we note that the 
analyses from Sections 4 and 5 carry through essentially unchanged for the model 
described above. 

The argument from Section 3, however, breaks down, because a particle moving 
through an environment of constant density m < 1 will have mass which does not 
decay exponentially, even when the loss of mass due to branching is taken into 
account. Instead, such a particle will (at large times) have a mass which is random 
and typically of order 0(1 — m). 

Because of this, the existing argument only establishes Proposition 3.1 in a 
highly weakened form, with the condition m > 1 rather than m > 0. (It is possible 
to do very slightly better, by considering a variable bound m = m{t). One can then 
take m(t) < 1 if 1 — m{t ) decays sufficiently quickly, but the pain-to-gain ratio in 
writing down such an argument in detail does not seem favourable.) But m € (0,1) 
is the really meaningful region. Proving a genuine analogue of Proposition 3.1 for 
this model seems to us the only missing step to a proof of Theorem 1.1 for the 
modified dynamics. 

• In the variant just described, one intriguing possibility is that there may now be 
particles with mass 0(1) at large times. If there are, they will be found near the 
front, since that is where they can find food. Do they really exist? 

• More generally, one may take 

Mi(t) = exp( [ {a - bC{s , Xj(s)))ds). 

Jo 

This looks, heuristically, like some sort of spatial logistic growth [ 8 , 15]. It may be 
interesting to investigate what different behaviours can occur as the parameters a 
and b are varied. 
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Appendix A. Estimates for the Upper Bound 

We first turn to the proof of Lemma 3.3. The proof relies on the following sample path 
estimate for Brownian motion. 
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Lemma A.l. Suppose f : [0,t] —» R and L : [0,f] —>• (0, oo) are twice continuously 
differentiable functions, with /(0) < 0, /(0) + L(0) > 0 and f increasing. We assume that 
there exists a constant Q > 0 such that 

\L'(0)\L(0) + \L'(u)\L(u) + [ U \L"(s)\L(s) ds + /“ |/"(s)|L(a) ds - |L'(0)|/(0) < Q 

Jo Jo 

/or all 0 < u < t, which we call Assumption (A). Then there is a constant M(Q) such that 
for 0 < p < 1. 

p { B ( s ) - f(s) € (0, L(s))Vs < u, B(u) - f(u) € (pL(u), L(u))} 

pu pu 2 

< M(Q) exp(—J / f'(s) 2 ds-j ^ ds - pf'(u)L(u) - f'(0)f(0) + \ log L(u)). 

Jo Jo 

This result is obtained by combining Proposition 4 and Lemma 7 in [20] to cover the 
two cases J Q U ds > 1 and J 0 “ -gJp ds < 1. In order to apply Lemma A.l, we exploit 
the existence of a solution to an integral equation; such a solution is used for a related 
purpose in Section 3.4 of [13]. 


Lemma A.2. For c < c*, there exists a constant a > 0 such that the equation 

-2 


l(s) = a + cs 1 / 3 — 


7r 


2y/2 


r s i 

Jo 


■ du 


(26) 


/ias a continuous solution on [0,1] which is twice continuously differentiable on (0,1) with 
l(s ) > 0 for all s € [0,1) and with 1(1) = 0. 


The lemma follows from Propositions 3.2 and 3.6(iii) of [13]. More precisely, in those 
lemmas there is a variance term a, and the value analogous to c* is a c = §(37r 2 cr 2 ) 1//3 . 
Taking o = l/\/2 yields the above formulation. 


Proof of Lemma 3.3. Fix t > 0 large. Since the integral on the RHS of (26) is non-negative, 
l(s) < c + a for all s € [0,1]. Since 1(0) = a > 0, we can fix (3 € (0, min(a 3 /8, a 3 /(8c 3 ), 1)) 
sufficiently small that l(s) > a/2 for s < /3. Let 


L(s) = f 1 / 3 



/ (s + /3t)u t \ 
V t + (3t ) 


(27) 


for 0 < s < t, where ut = infju € [0,1] : l(u) < 2 1 1 / 12 }. Note that ut —>• 1 as t —>• oo. 

We will prove that the lemma holds for above choice of (3 and with the function A(s) = 
L(s) — Kt 3 / 6 , provided K is sufficiently large. We must thus verify that A satisfies the 
requisite properties, and prove the bound (3). Write f(s) = b(s) — Kt 1 /® = \/2 .s — c(s + 
fit) 1 / 3 — Kt 1 /®; then to prove (3), it suffices to show that for u € [0, t] and all x € [0, L(u)), 

P (J?(s) € (f(s), f(s) + L(s)) Vs < u, B(u) > f(u) + x} < exp(-u— t 1 / 3 /K+\/2(L(u)—x)). 

(28) 

We establish this by applying Lemma A.l with the above functions / and L. We next 
derive the properties of /, L, L' and L" which we require to do so. 

First, note that /(0) = —c(f3t) 1 / 3 — Kt 1 /® and L(0) > t 1 / 3 l((/3ut)/( 1 + /?)) > at 1 / 3 /2 by 
our choice of /3, so since f3 < (a/ (2c)) 3 we have L(0) + /(0) > 0 for t sufficiently large. For t 
sufficiently large, / is increasing on [0, t ], and / and L are twice continuously differentiable 
(since l is C 2 on (0,1)). 

We assume t is sufficiently large that (1 + f3) 1 / 3 u t < 2. Then for s € [0, £], L(s) > 

2(1 + ( 0) 1 ' /3 'u/ 1 / 3 f 1 ' /4 > 2t 1 / 4 and L(s) < (1 + /3) 1 / 3 u/ 1 ^ 3 t 1,/3 (c + a) < 2(c + aft 1 / 3 , so for 
all s € [0, t ], 


21 1 / 4 < L(s) < 2(c + a)^/ 3 . 


(29) 
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Since l is C 2 on (0,1) and l(s) > 0 for s < 1, we can differentiate both sides of (26) for 
s € (0,1) to give 


l'(s) = -cs 2 / 3 — 

3 2y/2 Ks) 2 


TV 


1 


(30) 


Hence L is differentiable on [0, t] with 


L\s) = t 


- 2/3 


1 + P 

u t 


- 2/3 


, f u t {s + pt) 
t + pt 


(31) 


(32) 


Also, for 2^4 < u < ut, by (30) and the definition of ut, 

|/'HI < \c(pu t )- 2 /\ 1 + /3) 2/3 + < 2£ 1 / 6 , 

for t sufficiently large, so for all s € [0, t] we have 

|Z/(s)| < u t 2/3 (l + /?)— 2 / 3 t —2 / 3 2t 1 / 6 < 2£” 1/2 . 

This is a convenient moment to verify that the function A(s) = L(s) — ATf 1 / 6 has the 
requisite properties. By (29), for t sufficiently large, A : [0, t] —>• [t 1 / 4 , 2(c + a)t 4 / 3 ]. Also 
A(t) = L(t) — AT 1 / 6 < 2(1 + /3) 1 / 3 i^~ 1 / 3 t 1 / 4 < 4f 1 / 4 for t sufficiently large. Finally, by 
(32), |A'(s)| < 2t -1 / 2 < 1 for all s € [0,t], once again for t sufficiently large. 

Proceeding with the proof of (28), we now check that Assumption (A) holds for our 
choice of / and L, for some constant Q which does not depend on t. For t sufficiently large, 
by (29) and (32) we have sup sg [ 0if ] \L' (s)L(s)\ = 0(t _1 / 6 ), and also |L'(0)/(0)| = 0(t -1 / 6 ). 

By the definition of L in (27), for s € [0, t] we have 


|L"(s)|L(s) = t 


- 4/3 


1 + P 

u t 


- 4/3 


u t (s + pt) 
t T pt 


u t (s + pt) 
t T pt 


For < u < ut we also have |/'(u)| < 2t 4 / 6 by (31) and 2f -1 / 12 < l(u) < c + a by 
the definition of ut- By differentiating (30) we obtain that, uniformly over u in the above 
range, 

|i"(a)|((u) < |ar s/3 ((«) + • 

and hence sup s6[0)t ] \L"(s)\L(s) = 0(£ -1 )- Finally, sup se [ 0)t] /"(s) = sup sg[0it] §c(s + 
Pt)~ 5 / 3 = 0(t~ 5 / 3 ), which with (29) yields sup sg [ 0j t] \f"(s)\L(s) = 0(t~ 4 / 3 ). Thus, As¬ 
sumption (A) holds for some fixed constant Q not depending on t. 

Having verified the conditions of Lemma A.l, we now show that the bound from that 
lemma indeed implies (28). 

For 0 < u < t, f'{u) = y/2 — |c(s + pt)~ 2 / 3 so 

— \ f f\s) 2 ds = —u + V2c{u + pt) 1 / 3 — v / 2c(/3£) 1 ^ 3 + |c 2 (u + /It)” 1 / 3 — |c 2 (/3t)” 1 / 3 

Jo 

<-u + V2 cu 1/3 + 0(t” 1/3 ). 

Also by the definition of L in (27), 

1/3 


vr 2 r 

Jo 


L( s ) 


du = 


(l 


(u+P t)u t 
t+pt 


Put 

1 +p 


l(s)‘ 


■ ds 


We chose P sufficiently small that for s < ifjj, l(s) > a/2. Therefore 
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since we also chose /3 < a 3 / 8. It follows by (26) that 

V r ' in > f 1±^) 1/3 (V + „ ( < u ■+ ) 

2 7o L(s) 2 \ ut ) y2 \ t + fit ) 

> (y/2/2 )at 1/3 + V2cu l/3 - y/2L(u), 


1/3 


/ {u + f3t)u t \\ 

V t + pt )) 


where the second line follows by the definition of L in (27). Hence for u < t and p € [0,1], 

/*u /*u 2 

exp(-W f\s) 2 ds — ™ ds - pf'{u)L(u) - /'(0)/(0) + \ log L(u)) 

Jo Jo z ^{s) 

< exp(-u - (y / 2/2)at 1 / 3 + y/2(l - p)L(u) + V2cf3 1/3 t 1/3 + 0(t 1/6 )) 

< ( M{Q))~ l exp(— u — t 1 / 3 /K + \/2(l — p)L(u)). (33) 


The last inequality holds for all large t , provided K is sufficiently large that ( \/2/2)ot — 
\j2cfi 1 ! 3 > 1/K] this is possible by our choice of ft. 

Writing p = x/L(u), then 

P {B(s) € (f(s), f(s ) + L(s)) Vs < u, B(u ) > f(u) + x} 

=P {H(s) - f(s) € (0, L(s))Vs < u, B(u ) - f(u) € (pL(u), L(u))} 

< exp(— u — t 1 / 3 /K + V2(l — p)L(u )) 

for some iV > 0 by Lemma A.l and (33). This establishes (28) and completes the proof. □ 


We now turn to the proof of Lemma 3.5, during which we will use the following fact. 

Fact A.3. Let (W(u), 0 < u < 1) be either Brownian excursion or Brownian meander, 
and let N be a standard Gaussian. Then 

P |max W ( u ) > cc j> < 4P {N > x/4} . 

Proof of Fact A.3. Write B , B me , B ex and B hl for Brownian motion, meander, excursion, 
and bridge, all of length one. In what follows, maxima are always over u € [0,1] even if 
this is not explicitly written. 

We have 

maxB ex = max5 br — minH br < 2max|H br | = 2max|H(u) — uB{ 1)| < 4max |£>|, 

U<1 

so by the reflection principle, 

P {maxB ex > x} < P {max \B\ > x/4} < 2P {rnaxH > x/4} = 4P {N > x/4} . 

Next, let ( Z(u),u > 0) be a standard Bessel process. Since B me is Brownian motion 
conditioned to stay positive until time one, and Z is Brownian motion conditioned to stay 
positive for all time, it follows straightforwardly that max u <i Z(u) stochastically dominates 
max B me . By Pitman’s 2 M — B theorem (see [19]), we have 

ma xZ(u) = max | B(u) — 2 inf Bit)\ < 3max \B\, 

U<1 u <1 t<u 

and it follows as above that 


P {max5 me > x} < P {max \B\ > x/3} < 4P {N > x/3} . 


□ 


Proof of Lemma 3.5. First, for a standard Gaussian N we have P \N > cj < —l==^e c ~I 2 

L J v 27 T c 

for all c > 0. Using this bound, Fact A.3 and Brownian scaling, for each i we obtain 


max Xj(u) > x > < 4P 

U<ti 1 


N > 


x 


At. 


1/2 


< 


1 16/ 


1/2 


s —x 2 /16 ti 




< 


8 1 


1/2 


0 —x 2 /16 U 


X 
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Provided that x > 8T 1 / 2 , a union bound then yields 


P 


< max max X; (u) > x 

{ i> 1 u<ti 



e ~x 2 /16ti 


Finally, the function /(a) = e x2 ^ a is convex for a € [0, x 2 /2], and /(0) = 0, so if x 2 > 32 T 
then for each i, /(16L) < (tj/T)/(16T). Hence 


P 


< max max X;(u) > x 

( i> 1 «<i,; 


< e-* 2 / 16T 


□ 


Appendix B. Probability tail bounds 


We first state a Bernstein-type inequality due to Colin McDiarmid. 

Theorem B.l ([17], Theorem 2.7). Let Xi,... ,X n be independent with X^ — EX^ < 1 
for each k. Write S n = Y^k=i -^fc an d fix V > Var {<S n } = ]P/ =1 Var {-P*}- Then for any 

*^+**H^r <{*)'■ 

The first inequality is the heart of the theorem; the second is easy and is included to 
simplify an application of the theorem. The next lemma provides upper tail probability 
estimates for weighted geometric sums. 

Lemma B.2. Fix e < 1/2 and let ( Gi,i > 1) be i.i.d. Geom(l — e). For any n and any 
non-negative real numbers r\,...,r n with rriaxr,;/ ]P r* < 1/P, for all 6 > 0, 

{ n n \ 

Y r i°i > (1 + 5 ) n \ < 2(2 1+s e s ) v . 

i=l i=1 J 


Proof of Lemma B.2. Let Gj = Gj — 1 and pj = r 3 j ri. Then we must bound 

p x>4>4 , 


. i= 1 


under the assumption that ma XiPi < 1/P. First note that for c such that ee c < 1, 

1 - e 


E 


exp cP • y ^PiGj 


i= 1 


n 


1 - ee cV Pi 


For c > 0, the latter product is maximized (subject to the constraints that maxjpj < 1/P 
and that YliPi = 1) when pi = 1/|"P] for [P] values of i and Pi = 0 otherwise. We thus 
obtain 

(1-0™ 


E 


exp(cP • y^PiGj) 


i=l 


< 


(1 - ee c )TP 


For any non-negative random variable, P {X > <5} < e c<5 ' Ee cVX ; taking e c = (2e) 1 
yields 


p &g,>4< 


-cSV 


(1-C 


\V 


. i —1 


(1 - ee c ) v+l 


= 2 


(2(l-e)(2e 


v 


< 2(2 i+ V 


l+<5 J\V 


□ 
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